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We consider the Rabi Hamiltonian which exhibits a quantum phase transition (QPT) despite consisting only 
of a single-mode cavity field and a two-level atom. We prove QPT by deriving an exact solution in the limit 
where the atomic transition frequency in unit of the cavity frequency tends to infinity. The effect of a finite tran¬ 
sition frequency is studied by analytically calculating finite-frequency scaling exponents as well as performing 
a numerically exact diagonalization. Going beyond this equilibrium QPT setting, we prove that the dynamics 
under slow quenches in the vicinity of the critical point is universal, that is, the dynamics is completely charac¬ 
terized by critical exponents. Our analysis demonstrates that the Kibble-Zurek mechanism can precisely predict 
the universal scaling of residual energy for a model without spatial degrees of freedom. Moreover, we find that 
the onset of the universal dynamics can be observed even with a finite transition frequency. 


Introduction.- Universality plays a key role for our un¬ 
derstanding of quantum phase transitions (QPT) in interact¬ 
ing quantum systems ID. While the concept of universal¬ 
ity is well-established in equilibrium QPT, the question to 
what extent the concept of universality could be extended to 
non-equilibrium dynamics of QPT remains largely to be ex¬ 
plored ei a. For a slow quench across a QPT, the closing 
spectral gap at a critical point leads to a breakdown of the adi- 
abacity regardless of the quench rate. The scaling of defect 
formation has been shown to be entirely controlled by criti¬ 
cal exponents and quench rate through a successful applica¬ 
tion of the Kibble-Zurek mechanism (KZM) mm, originally 
developed for classical phase transitions, to QPT in short- 
range interaction models f8 UT2l . However, whether this scal¬ 
ing holds for fully-connected models ED, which lack spatial 
degrees of freedom, such as Dicke m or Lipkin-Meshkov- 
Glick (LMG) model (SI remains an open problem CMS). 

The Dicke model considers a system of a quantized single¬ 
mode cavity field uniformly coupled to N two-level atoms. 
It exhibits a superradiant QPT in the thermodynamic limit 
(N -A oo) QSH22. While tremendous efforts have been 
devoted to understand the QPT of the Dicke model both 
in and out of equilibrium Ifl9ll27l . a criticality of the Rabi 
model II27U331I . the most simplified version of Dicke model 
with N = 1 , has been hitherto largely overlooked. Having 
only two constituent particles, the Rabi model is far from 
being in the thermodynamic limit where a QPT typically 
occurs; however, a ratio of the atomic transition frequency 
Q to the cavity field frequency u >0 that approaches infinity, 
Q/wo -A 00, can play the role of a thermodynamic limit f27t 
that allows the spectral gap to be precisely closed at the criti¬ 
cal point ID. 

In this letter, we firstly establish the theory of equilibrium 
QPT of the Rabi model. At the core of our analysis is a low- 
energy effective Hamiltonian that is valid for Q/wo 1 and 
becomes exact in the Q/wo —> 00 limit. We derive an ex¬ 
act solution for eigenstates, an energy spectrum, expectation 
values of relevant observables, and critical exponents in the 
Q/w0 -A 00 limit. Our solution shows that there exists a crit¬ 
ical atom-cavity coupling strength g c beyond which the /Q 
parity symmetry is broken and the cavity field is macroscop- 


ically occupied. Further, the effect of a finite value of fi/tu0 
on the QPT is analyzed in the spirit of the finite-size scaling 
analysis. The leading order corrections to the ground state en¬ 
ergy, the excitation energy, the average photon number, and 
the variance of cavity field quadratures at the critical point 
are derived analytically, from which finite-frequency scaling 
exponents are obtained. We also perform an exact diagonal¬ 
ization and find an excellent agreement with analytical results. 

Our establishment of the equilibrium QPT allows us to in¬ 
vestigate the universality in the dynamics of the Rabi model. 
Particularly, we are interested in quench dynamics where the 
system is initially prepared in the ground state and the control 
parameter g is tuned towards to the critical point linearly in 
time with a quench time r q starting from g = 0 mmm 
[36l . On the one hand, we solve the dynamics exactly in the 
Q/Wi -A 00 limit, and calculate the residual energy as a mea¬ 
sure of the degree of non-adiabacity, which shows a power- 
law scaling with the quench time r q . On the other hand, we 
obtain such a scaling solely from the critical exponents found 
in the first part of the letter. To this end, we apply KZM to 
the adiabatic perturbation theory I9l f34)436ll and the dynamical 
critical function method EO , independently. Both approaches 
give rise to the same universal scaling that precisely predicts 
the exact dynamics, demonstrating that the KZM can lead to 
a universal dynamics for a model without spatial degrees of 
freedom. 

Finally, we consider the same quench dynamics with a fi¬ 
nite value of n/u 0, and show that, as one decreases the ratio 
Q/wo, there is a crossover from the universal scaling to the 
r“ 2 scaling, a typical scaling of the adiabatic dynamics with 
a finite quench time for a gapped system iflfil [34] l35l. We 
identify a range of quench times which leads to dynamics that 
closely follows the universal scaling, and show that the onset 
of the universal dynamics can be observed for a finite Vt/uj o- 
The crossover from the universal to the r~ 2 scaling is also 
observed in the fl/co 0 -A 00 limit by ending the quench of 
the control parameter g below the critical point. It demon¬ 
strates that the spectral gap opening due to finite Sl/w0 has the 
same effect as ending the quench below the critical point in 
the Q/wo -A 00 limit. 

Quantum phase transition.- We consider the Rabi Hamil- 
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tonian l37l 

H Rabi = uqo 1 a + j&z ~ A(a + a t )cr x (1) 

where a x>z are Pauli matrices for a two-level atom and a (o') 
is an annihilation (creation) operator for a cavity field. The 
cavity field frequency is ojq, the transition frequency Cl, and 
the coupling strength A. We denote (4-)) as eigenstates 
of o z , and |to) the eigenstate of a' a. The parity operator, 
II = e l ' K ( a a +j( 1 + cr *)) ) which measures an even-odd parity of 
total excitation number, commutes with // Ka bi- The Z) par¬ 
ity symmetry has been shown to be sufficient for the model 
to be integrable 1351 : however, a lack of a closed form solu¬ 
tion makes the approach in Ref. Il38l not directly applicable to 
investigate the QPT. 

In the CI/ojo —>• oo limit, we firstly find a unitary transfor¬ 
mation, U = exp[^(a + a/)(cr+ — rr_)], which makes the 
transformed Hamiltonian, V' H Rab iU, free of coupling terms 
between spin subspaces TQ and 'H~ . Upon a projection onto 
T~L^, i.e., H np = (J, \WH Rab iU\ 4-), we obtain an effective 
low-energy Hamiltonian, 

7 t t ^0 9 t , -j-\2 ^ 

Hnp — U 0 CL CL “ \CL -\- CL ) —j (2) 

where g = 2A/ \Joj 0 il |39l . Eq. (J2ji can be diagonalized to 
give H np = Cnptfb — Cl/2 with e np = wo\/l— g 2 , which is 
real only for g < 1 and vanishes at g = 1, locating the QPT. 
The low-energy eigenstates of // Kah j for g < iare 1 4>™ p {g)) = 
\ m ) I4-) with 5[a:] = exp[f (a t2 -a 2 )] and r np (g) = 

-jbi(l -9 2 )- 

The failure of Eq. Q for g > 1 suggests that the number of 
photons occupied in the cavity field becomes proportional to 
CI/ujo so that the higher order terms cannot be neglected, i.e., 
superradiance occurs; it also suggests that V\, is no longer the 
low-energy subspace. In order to properly capture the low- 
energy physics, we transform // Rah , of Eq. 0 by displacing 
the cavity field a, i.e., JT Rab i(±a ff ) = [±a g ]H Rabi D[±a g \ 

with T>[a\ = e Q ( a+-a ) and a g = \J (ff 4 — 1), which 
reads 

H Rabi (±a g ) = w 0 o t a + - A (a + a))r^ + uj 0 a 2 g (3) 

where r± = If*) (^H^) 4*1 = Eq. (g 

has the same stmcture as Eq. 0 with rescaled frequencies 
A = 4/°° and Cl = g 2 CL Therefore, by employing the same 
procedure used to derive H np , we find an effective Hamilto¬ 
nian of the Rabi Hamiltonian for g > 1 from Eq. 0, 

H sp = w 0 a)a - ^(a + a t ) 2 - j(g 2 + g~ 2 ), (4) 

whose excitation energy is found to be e sp = uiq \/l — g~ 4 , 
which is real for g > 1. Note that two independent choices 
of a = ±cr 9 in Eq. ([3]) lead to an identical spectrum. The 
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FIG. 1. Top panel: exact solutions of the Rabi model in the Cl/uj 0 —> 
00 limit as a function of the dimensionless coupling strength g/g c for 
(a) the rescaled ground state energy eG (solid) and d 2 ec/d 2 g (red- 
dashed), (b) the excitation energy e (solid) and the energy difference 
between the ground and the first excited state (red-dashed) showing 
the ground state degeneracy for g/g c > 1, and (c) the variance of 
position Aa; (solid) and momentum Ap (red-dashed) quadrature of 
the cavity field, and AxAp (dotted). In (b) and (c), the scaling rela¬ 
tion near the critical point is indicated. Bottom panel: a leading order 
correction for finite fi/wo at g = g c for Ap, e, the order parameter 
n c , and ec from top to bottom, respectively. The analytical results 
(lines) predict precisely the exact diagonalization results (points) for 
all observables. The finite-frequency scaling exponents for each ob¬ 
servable are indicated. 


low-energy eigenstates of // Ra b, for g > 1, WM)± = 

V[±a g ]S[r sp (g)} \m) 14*} where r sp (g) = -\ ln(l - g ~ 4 ), 
are, therefore, degenerate; they also have a spontaneously bro¬ 
ken parity symmetry, as evident from the non-zero coherence 
of the field (a) = ±a 9 . The higher order corrections in Eq. (|2j) 
and (|4j» vanish exactly in the Q/ojo —» oc limit. Therefore, 
H np and H sp are the exact low-energy effective Hamiltonian 
for the normal phase (g 1) and superradiant phase (g > 1), 
respectively, for which the subscripts np and sp stand. See 
Ref. ROl for a detailed derivation of the effective Hamiltonian 
and its solution. 

Our exact solution shows that the superradiant QPT occurs 
at the critical point g c = 1. The rescaled cavity photon num¬ 
ber n c = ^ (a^a) is zero for g < g c and n c = {g 4 - g/,)/^g 2 
for g > g c \ thus, n c is an order parameter. The rescaled 
ground state energy, e G (g) = E G (g ) is -w 0 /2 for g < g c 

and —cuo(g 2 + g~ 2 )/A for g > g c . While e G (g) is continu¬ 
ous, d 2 e G (g)/d 2 g is discontinuous at g = g c , revealing the 
second-order nature of the QPT [Fig. |T| (a)]. Near the criti¬ 
cal point, the excitation energy in both phases, e np and e sp , 
vanishes as e(g) oc \g — g c \ zv with zv = 1/2 [Fig. 0 (b)], 
where v (z) is the (dynamical) critical exponent. Meanwhile, 
the variance of position quadrature of the field x = a + 
diverges as Ax(g) oc \g — g c | -1 / 4 oc e~ 1//2 , from which we 
find that z = 2 and u = 1/4 [Fig. [T] (c)]. While we have de¬ 
fined the critical exponents z and v separately by noticing that 
Ax plays an analogous role of the diverging length scale in 
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extended quantum systems ID, only is the product zv an im¬ 
portant exponent in the following analysis. The critical point 
also accompanies an infinite amount of squeezing in the mo¬ 
mentum quadrature p = i(a' — a) so that it remains in the 
minimum uncertainty state for any g, i.e., Ax(g)Ap(g) = 1 
[Fig. |T](c)]. 

Finite-frequency scaling - We complete our study of the 
equilibrium QPT by investigating the, finite-frequency effect. 
Firstly, we derive a leading order correction to the exact effec¬ 
tive Hamiltonian. To this end, we find a unitary transformation 
U n = exp[(^(a + a t ) - |^(a + a t ) 3 )(cr + - cr_)] ofEq. |IJ 
that decouples the 77, and 'H- subspaces up to fourth order in 
X/fl and project to Ti\. to obtain BOl 


fjQ _ IT 

11 np — 12 np 


n 4 t i 2 n 2 t t 2 

+ ^0( a + a t)4, 9_ff0 
160 1 j + 4H ’ 


(5) 


where the leading order correction adds a quartic potential for 
the cavity field. Although H^ p is not exactly solvable, a vari¬ 
ational method can be used to derive analytical expectation 
values |40| . We find that, at the critical point, the excitation 
energy vanishes and the characteristic length scale diverges 
with a power-law scaling. 


e 9c (Q/cj 0 ) = ui 0 


/2H\ 1/3 
\3cjq/ 


AXg c (fl/UJ 0 ) 


/ 2fl \ 1/6 
\3wo/ 

( 6 ) 


In addition, the leading order correction for ec and n c 
are given by eG, go {Xl/co o) = (w 0 /4)(2fi/3u;o) _4//3 and 
n Ct g c {fl/uif) = l/6(2fl/3wo) _2//3 . The exponents of these 
scaling relations, the finite-frequency scaling exponents, are 
found to be the same as the finite-size scaling exponents of 
corresponding observable for the Dicke model RTl and LMG 
model (4211431 . which also have the same critical exponent z 
and v (4411451 . We perform an exact diagonalization of Eq. (0 
and show that the numerically obtained scaling exponents pre¬ 
cisely match the analytical results [Fig. |T](d)] . 

Universal scaling for adiabatic dynamics.- Having estab¬ 
lished the equilibrium QPT of the model, we are now able 
to investigate the dynamics of the QPT. We consider a proto¬ 
col where the control parameter g is changed linearly in time, 
g(t) = gft/rq , with gf being the final value. The system is 
initially in the ground state. As gif) approaches the critical 
point, the vanishing spectral gap makes the relaxation time 
of the system diverge, inevitably creating quasiparticle exci¬ 
tations irrespective of how large the quench time r q is. Ap¬ 
plying kzm EEHEQ, we define a time instant t that divides 
the dynamics into the adiabatic and impulsive regime from 
g 2 (t) = f]{f), where the accessible energy gap g is given 
as ?/ = 2e np for g < g c due to the parity symmetry. From 

e np = woy/ 1 - 9 2 - we find 9 ~ 9c ~(±s/2ua gO) 
where the coupling instant g = gif = t) moves away from the 
critical point as one decreases the quench time so that the im¬ 
pulsive regime widens. Note that we only consider gif) < g c 
for a simplicity (34U361 . 
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FIG. 2. Residual energy E r as a function of the quench time r q ob¬ 
tained by solving a nearly adiabatic dynamics for (a) different values 
of the final coupling strength g f ranging from gj = Qhg c to gf = g c 
(from bottom to top) in the Q/ujo —> oo limit and (b) different ratios 
H/wo ranging from fl/tuo = 10 2 to Q/luo —» oo (from bottom to 
top) with a fixed final coupling strength gf = g c . For gf = g c in the 
0,/ujo —» oo limit, it precisely follows the universal scaling relation 
(solid line), predicted by the Kibble-Zurek mechanism. Both moving 
gf away from g c and reducing the ratio fl/uto result in a crossover 
from the universal scaling to rf 2 scaling (dashed line). 


The wave function at time t can be expressed in terms 
of the instantaneous eigenstates of H np (g(f)), i.e., l’P(i)) = 
^ m c m (t) < S[r Tl p(f)] |to). Then, we apply the adiabatic per¬ 
turbation theory (9l l34ll36l to calculate the residual energy E r 
at the end of the quench, which measures the degree of non- 
adiabacity, defined as E r = (4'(T 9 )|7T Tlp (()(r <? ))|'I'(r 9 )) — 
Ec(g(Tq)). For a protocol that stays in the adiabatic regime, 
i.e., gf <C g, we obtain a scaling relation, E r ex r~ 2 l40l . 
which is a typical scaling for the adiabatic dynamics with a 
finite quench time for a gapped Hamiltonian. If the protocol 
involves the impulsive regime, gf ~ g, we find that the resid¬ 
ual energy follows a universal scaling relation, 

E r oc T f zv ^ zv+1 \ (7) 

that is, E r oc since zv = 1/2 (40l . A different way to 

predict the universal scaling of E r based on KZM is to use the 
dynamical scaling function approach 03, which expresses 
the scaling relation in terms of the finite-frequency scaling ex¬ 
ponents. We confirm that it predicts the same universal scaling 
relation as in Eq. ([7]) l40l . 

For short-range interaction models, the residual energy due 
to a slow quench stems from spatial defects in order param¬ 
eter across a QPT, whose scaling has been successfully pre¬ 
dicted by KZM EUSE]. However, it is not clear whether 
KZM can predict the scaling of the residual energy in fully- 
connected models due to their lack of spatial degrees of free¬ 
dom. In fact, although the same scaling relation with Eq. (j7l) 
has also been predicted for the Dicke and LMG model [ 171 . 
a numerical calculation with a finite-size LMG model shows 
a significant discrepancy with the universal scaling as it esti¬ 
mates E r ex T q 3 ^" El , raising a doubt on the applicability of 
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FIG. 3. Exponents of power-law scaling rf for the residual energy 
E r presented in Fig. [2] (b). (a) Fits obtained for different ranges of 
the quench time [r 9 i, r q 2 \- The values for r ql ( q2 ) are indicated in the 
figure. The exponent g converges to the universal scaling exponent 
— 1/3 for finite £1/ujq. (b) Fits obtained for a range of the quench 
time [r 9 Ar 9 , r q /Ar q ] as a function of r q with a fixed log-scale in¬ 
terval log 10 Ar, = 6.25 x 10 -2 . The crossover from g = —2 to 
g = — 1/3 as one increases G/cjo is clearly demonstrated. 


the KZM to the fully-connected models 03 . Strictly speak¬ 
ing, one has to solve the dynamics exactly in the thermody¬ 
namic limit for the LMG or Dicke model, or equivalently in 
the G/wo —> oo limit for the Rabi model to test the validity 
of the universal scaling relation, which is accomplished in the 
following section. 

Exact solution for adiabatic dynamics - The exact low- 
energy effective Hamiltonian given in Eq. 0 allows one 
to numerically solve the slow quench dynamics of the Rabi 
model, which involves only a small number of quasiparticle 
excitations, in the G/wo —> oo limit. The equation of mo¬ 
tion is given as idfj(t) = [a H (t), H nPiH (t)\, where the sub¬ 
script H indicates the operators in the Heisenberg picture. 
We express the cavity field operator at time t as anit) = 
u{t)a( 0) + t)*(f)a^(0) with an initial condition w(0) = 1 and 
n(0) = 0, and \u(t)\ 2 — |i;(t)| 2 = 1, and derive coupled dif¬ 
ferential equations for u{t) and v (f), 


i du(t) 
wo dt 
i dv{t) 
wq dt 



u(t ) 
v(t) 


g 2 (t)_ 


it), 


g 2 (t ) 


it)- 


( 8 ) 


The residual energy in terms of u(t ) and v(t) is given by 


E r = u 0 \v(t)\ 2 


w 0 gj 
4 


\u(t)+v(t) | 2 


^np(gf) Wo 
2 


(9) 


In Fig. [2] (a), we plot E r at the end of the quench as a func¬ 
tion of r q for different values of the final coupling strength 
gf in the G/wq oo limit. For a protocol that ends right at 
the critical point, gj = g c , it precisely follows the universal 
scaling given in Eq. 0 It confirms that the nearly adiabatic 
dynamics of the QPT in the Rabi model can be completely 
characterized by the critical exponents alone, thus is univer¬ 
sal. We note that the saturation of E r observed for a short 
quench time, r q < 1/wo, corresponds to sudden quench dy¬ 
namics. As we change gj progressively away from the critical 


point, gy < g c , the universal scaling breaks down and the t~ 2 
scaling emerges, which is precisely the scaling predicted using 
the adiabatic perturbation theory in the adiabatic regime iia. 

We find that the leading order correction to the equation 
of motion for finite G/wo adds an additional term, f(u, v) = 
(3wo/4G)g 4 (t)(w + u)|u + v\ 2 , to the right hand side of both 
equations in Eq. (|8]i 1 401 . For a quench that ends at the critical 
point gf = g c , the leading order correction to the residual en¬ 
ergy adds an additional term, h(u,v ) = (3wo<? 4 )/(16G)|m + 
v\ A - ^(2G/3w 0 )" 1/3 , to Eq. ^ gO). In Fig.||(b), where 
we reduce the ratio G/wo from infinity to 10 2 for gf = g c , we 
observe a crossover behavior for the residual energy virtually 
identical to Fig. [5] (a). This is because a finite value of G/wo 
opens up an energy gap at gf = g c whose effect is equivalent 
to ending the protocol away from the critical point. 

An interesting aspect of the crossover behavior for the scal¬ 
ing of E r shown in Fig.[2](b) is that there is a range of quench 
time Tq at around r q £ [10,10 3 ] where the E r closely fol¬ 
lows the universal power-law even for finite values of G/wo- 
By closer inspection, we find fits for the slope of curves in 
Fig. [2] (b), which corresponds to the exponents of power-law 
scaling of E r , for a wide range of quench times. As shown 
in Fig. [3] (a), the exponents converge to the universal scaling 
exponent —1/3 as one increases the ratio G/wo, showing that 
the onset of the universal dynamics can be observed with finite 
G/cl>o- The convergence to the universal scaling implies that 
the energy gap whose scaling is given in Eq. (|6| is sufficiently 
small to drive the system into the impulsive regime so that the 
dynamics is strongly influenced by the nature of the critical 
point. As the energy gap widens for smaller values of G/wo, 
the influence of the critical points gradually vanishes, leading 
to a crossover to r“ 2 scaling. In Fig.pl(b), the crossover of 

_i/o 

the scaling from r q ' to r“ 2 is further elucidated by find¬ 
ing fits for much shorter interval of r q , which approximates 
the slope of the tangent line of graphs in Fig. |3](b). Varying 
gf in the G/wo —> oo limit shows identical features shown in 
Fig.0j4Ol. 

Conclusion - We have found an effective low-energy de¬ 
scription of the Rabi model that unveils the universality of the 
model both in and out of equilibrium. Our analysis shows 
that the superradiant QPT which has been primarily studied 
for systems of thermodynamically many atoms can as well 
be investigated with systems of a single atom. An important 
advantage of the reduced degrees of freedom is that solving 
the critical dynamics is more tractable; indeed, we have been 
able to report a first confirmation of the KZM prediction for 
a model without spatial degrees of freedom. Together with 
an impressive ongoing progress of technologies to realize the 
interaction between a two-level system and a single harmonic 
oscillator, we expect that the Rabi model can serve as an excel¬ 
lent platform to study equilibrium and non-equilibrium criti¬ 
cal phenomena. 
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